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ABSTRACT 

We show that the globular cluster luminosity function (GCLF) of the inner halo of 
the Milky Way is statistically different from the GCLF of the outer halo. We also find 
a similar difference between the inner and outer halo population of M31. We assert 
that this difference is evidence for some form of dynamical evolution of the cluster 
population and/or a dependence of GCLF shape on the environment in which the 
cluster population formed. We also find that the turnover luminosity of the GCLF is 
unaffected by these differences and further assert that this stability of the turnover 
luminosity affirms its usefulness as an indicator of cosmic distance. 
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1 INTRODUCTION 

The globular cluster system of the Milky Way has long been 
known to be a multiple component system that can be di- 
vided into disk and halo subsystems; in addition, the halo 
subsy stem can be further d ivided into an inner and outer 
halo (Searle and Zinn 197£). This division of the halo into 
inner and outer components has allowed the discrimination 
of the clusters along age boundaries on the basis of hor- 
izontal branch type (Lee et. al 1994), and the picture of 



inside-out galaxy formation now seems secure. For an excel- 
lent review of these issues see Zinn (1996) and other works 
in the same volume. 

Consistent with the emerging acceptance of an inside- 
out picture of galaxy formation has come the view that the 
current population of globulars is not the primordial distri- 
bution but r ather the surviving component of a once grande r 
population ( Fall and Rees 1977 ; Murali and Weinberg 1996 ). 
Most recently, a very appealing picture has been presented 
in which outer halo clusters do not suffer significant dis- 
r uption while those which a re interior to Rgc ~ 8kpc do. 
(Murali and Weinberg 1996). 

Given both that the outer halo clusters are younger and 
that their orbits make them less likely to have undergone ex- 
ternally imposed dynamical evolution, the outer halo cluster 
population should be most like that of a primordial system 
while that of the inner halo should be evolved. This reason- 
able view can be tested by a comparison of the globular clus- 
ter systems (hereafter GCSs) of the inner and outer halo. To 



allow an extension of this comparison to external galaxies, 
the comparison should be between observables which can be 
studied in remote systems. The globular cluster luminosity 
function (GCLF) is the obvious choice. 



2 THE MILKY WAY GCLF, INSIDE AND OUT 

Figure |l| shows (a) the GCLF for all of the halo clusters 
of the Milky Way, and the GCLFs for the inner (b) and 
outer (c) halo subsystems. (These data were taken from the 
McMaster University globular cluster database maintained 
by W.E. Harris^]). There are approximately 50 clusters in 
each of the subdivisions and so the differences between these 
two distributions are decidedly real. A Kolmogorov-Smirnov 
comparison of the two populations gives a probability of 
less than 5% that the two distributions come from the same 
parent population. This evidence strongly suggests that the 
inner and outer halo populations have different formation 
histories and/or different evolutionary histories. 
Fits of the function 



vlexp(-(m - M T o) 2 /(2cr 2 )) 
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to the three distributions are overlaid on the histograms 
in Figure with the numerical results shown in Table |l| 
These fits were performed using the NGAUSSFIT routine in 
STSDAS under the assumption of Poisson sampling errors. 
It is of interest to note that although the inner and outer halo 
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Table 1. Fits to the GCLFs for various components of the 
Milky Way GCS. These fits were performed using the NGAUSS- 
FIT package in STSDAS. The data are from the McMaster cata- 
logue. Note that the peak or turnover magnitude remains constant 
(within uncertainties) for all components in the Galaxy but the 
outer halo population is significantly broader. 
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Figure 1. (a) The luminosity function for the halo globular clus- 
ters of the Galaxy, (b) The luminosity function for the inner halo 
clusters, (c) The luminosity function for the outer halo clusters. A 
Kolmogorov-Smirnov test shows that there is a 3.7% probability 
that the two distributions are drawn from the same parent pop- 
ulation (ie. the hypothesis that these two distributions are from 
the same parent population is rejected). 



GCLFs have markedly different distributions, they have the 
same peak luminosity. This res ult appears to be in conflic t 
with recent dynamical models ( Murali and Weinberg 1996 ) 
which suggest that the peak value of the GCLF will become 
significantly fainter as the cluster population evolves, and 
that such evolution should be strongest for the inner halo 
populatieHr 



3 THE GCS OF M31 

The recognition of the two-component nature of the Milky 
Way halo cluster luminosity function immediately suggests 
an examination of the M31 GCLF under the same condi- 
tions. To allow such a comparison it is necessary to divide 



the M31 cluster population into disk and halo subsystems. 
The question of which M31 clusters are halo members and 
which are disk members is confused by the lack of spectro- 
scopic metallicities for a vast number of cluster candidates. 
Ideally a cut of metallicity similar to that used for the Milky 
Way would isolate the disk and halo subsystems, but in prac- 
tice such a straightforward approach is not possible. 

Previous authors have selected clusters as halo members 



based on their position on the sky (Reed et. al 1994). This 
selection criterion has the disadvantage of preferentially ex- 
cluding many inner halo clusters. To avoid excluding objects 
on the basis of position we select as our "halo sample" all 
clusters for which (B-V) < 0.8 in the Battistini et al. (1987) 
survey. This selection seams justified in that for the Galaxy 
there are no disk clusters bluer than (B — V) = 0.8 and so 
this cut should eliminate the majority of disk clusters from 
the M31 sample. Unfortunately this selection criterion also 
removes many true halo members from the sample; however, 
we find that there is no significant difference in parameters 
between a Milky Way halo sample selected as [Fe/H] < —0.8 
and one selected using (B — V) < 0.8. This culling also has 
the advantage of excluding background galaxies from th e 
Battistini et al. (1987) A,B sample (see ( Reed et. al 1992 ) 
Figure 8). 

In panel (a) of Figure we present the LF for all the 
globular cluster candidates which have (B-V) < 0.8 (from 
Battistini (1987)). In panel (b), we present the luminosity 
function for the subset of objects which lie within a projected 
radius of lOkpc (140 ar cmin for a true distance modulus of 



24.45 ( [lacoby et al. 1992) )) of the center of M31; in panel 
(c) we show the luminosity function for the objects beyond 
lOkpc. This figure qualitatively reveals the same separation 
of the GCLF into peaked and flat components, when selected 
on the basis of radius, as was seen for the Galaxy. We do 
not, however, consider this as strong a test as that afforded 
by the MW sample because of the uncertain extent to which 
incompleteness affects the luminosity functions in the sub- 
samples at the faint end. In addition, of course, there will be 
some contamination of the "inner halo" subsample by outer 
halo clusters projected into the 10 kpc circle. Also, the de- 
termination of halo/disk membership for individual clusters 
is uncertain and it is probable that there is contamination 
of the inner halo sample by disk clusters. 

Table |2| presents the results of fits of equation jj] to 
the various cluster populations in M31. As was found for 
the Milky Way GCS the outer halo clusters of M31 have 
a broader luminosity function, and once again there is no 
evidence for difference between the turnover luminosity of 
the inner and outer halo. A K-S comparison of the data sets 
gives a probability of 18% that the two samples are drawn 
from the same parent population. This provides marginal 
evidence for a multi ple component halo and is consistent 



with previous work (Ashman and Bird 1993 



Huchra et al. 



1991) which has shown that the distribution of halo clusters 



in M31 contains substructure. 

As for the Milky Way, the observed difference between 
the two M31 halo populations cannot be the result of dynam- 
ical evolution of the type predicted by Murali and Weinberg 
(1996) since their model implies that the inner population 
should appear fainter than that of the outer. There is no 
evidence for such a shift in turnover values. 
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Table 2. Fits to the GCLFs for various components of the 
Andromeda galaxy (M31). These fits were performed using the 
NGAUSSFIT package in STSDAS. The data are from Battistini 
et. al (1987). Once again, the outer halo population is somewhat 
broader than that of the inner halo sub-system. 



M31-GCS 
all halo clusters 
inner halo clusters 
outer halo clusters 



M TO 
17.45 ± 0.08 
17.30 ±0.10 
17.29 ±0.17 



VGCLF N 

0.76 ±0.05 161 

0.57 ±0.10 64 

0.94 ±0.13 97 
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ther assurance that GCLFs can be used as standard candles 
in cosmic distance determinations. The dependence of the 
shape of the GCLF on galactocentric radius suggests that 
future comparisons should be made between cluster popu- 
lations that only include clusters of comparable radii from 
their host galaxy centers. 

If dynamical effects are not found to be responsible for 
the dependence of GCLF shape on radius, then the envi- 
ronment in which clusters form is likely to be the deciding 
factor. 

Many thanks to S.T. Butterworth and D. Wing for their 
clarifying discussion of this issue. 
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Figure 2. (a) The luminosity function for the halo globular clus- 
ters of the M31. (b) The GCLF for the inner halo clusters, (c) The 
GCLF for the outer halo clusters. A Kolmogorov-Smirnov test 
shows that there is a 18% probability that the two distributions 
are drawn from the same parent population (ie. the hypothesis 
that these two distributions are from the same parent population 
is not supported). 



4 CONCLUSIONS 

Both within the Milky Way and, to a lesser extent, M31 the 
outer halo clusters are clearly more indicative of a broad lu- 
minosity function. This population may be, then, indicative 
of an initial broad mass function for globular clusters. 

The inner halo clusters demonstrate a remarkably 
peaked distribution. The difference between the inner and 
outer luminosity functions may indicate the fate that dy- 
namical evolution has in store for clusters formed near the 
centers of galaxies. 

The turnover luminosities of the inner and the outer 
cluster populations appear to be consistent. If this is the 
case, then dynamical models of cluster evolution will need to 
account for the preferential stripping away of clusters fainter 
and brighter than the turnover luminosity. This decoupling 
of mean luminosity from position in the galaxies gives fur- 
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